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(*~^ ' We examine weak-coupling theory for unconventional superconducting states of cubic or tetrahe- 

\ dral symmetry for arbitrary order parameters and Fermi surfaces and identify the stable states in 

£SJ ■ zero applied field. We further examine the possibility of having multiple superconducting transi- 

tions arising from the weak breaking of a higher symmetry group to cubic or tetrahedral symmetry. 
Specifically, we consider two higher symmetry groups. The first is a weak crystal field theory in 
which the spin-singlet Cooper pairs have an approximate spherical symmetry. The second is a weak 
spin orbit coupling theory for which spin-triplet Cooper pairs have a cubic orbital symmetry and an 
approximate spherical spin rotational symmetry. In hexagonal UPt3, these theories easily give rise 
to multiple transitions. However, we find that for cubic materials, there is only one case in which 
two superconducting transitions occur within weak coupling theory. This sequence of transitions 
does not agree with the observed properties of PrOs4Sbi2. Consequently, we find that to explain 
two transitions in PrOs4Sbi2 using approximate higher symmetry groups requires a strong coupling 
£j . theory. In view of this, we finally consider a weak coupling theory for which two singlet represen- 

tations have accidentally nearly degenerate transition temperatures (not due to any approximate 
symmetries). We provide an example of such a theory that agrees with the observed properties of 
PrOs 4 Sbi 2 . 
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It has been observed that a wide variety of heavy fermion superconductors appear to undergo multiple phase 
transitions within the superconducting state. These materials include UPt^i, Ui-^Th^Beia^ and PrC^Sbisf^- 
The fact that multiple transitions are being observed in such a significant fraction of the total number of heavy 
fermion superconductors discovered is presumably providing a valuable insight into the nature of the superconducting 
l— ~~ '■ state. Among the heavy fermion superconductors showing multiple phase transitions, UPt 3 has been most extensively 
I ' studied. It has a hexagonal point group symmetry and shows two transitions with a small (about 10%) splitting 
, of the superconducting transition temperature. Three main approaches have been used to theoretically model the 
phase diagram of UPt3 and it would be of interest to see if any of these approaches apply to the other materials. 
The first approach uses a weak symmetry breaking field (SBF) to break the hexagonal symmetry. This field lifts the 
degeneracy of a multi-component order parameter leading to two transitions^. The origin of this symmetry breaking 
field has been questioned^ and other proposals have emerged. Zhitomirsky and Ueda have postulated a weak crystal 
field model in which spin singlet Cooper pairs experience an approximate spherical symmetry which is weakly broken 
to hexagonal symmetry^. Along similar lines, Machida and co-workers have postulated a weak spin orbit coupling 
for spin-triplet Cooper pairs^ . The third approach uses a phenomenological model which considers two different 
irreducible representations of the hexagonal point group that accidentally have nearly the same T,-^!^ . In all these 
cases, it is possible to develop microscopic theories based on weak coupling theory that give rise to the two transitions 
I ■ and such theories have been useful in developing a qualitative understanding of this superconductor—. 

Among the materials with cubic or tetrahedral point group symmetry: Ui_ x Th x Bei3 and PrOs4Sbi2 have shown 
multiple superconducting transitions. Both these materials have been the subject of phenomenological studies. 
Ui_ x Th x Bei3 has been studied using a phenomenological approach by Sigrist and Rice 13 while PrOs4Sbi2 has been 
examined phenomenologically by Goryoi*, Ichioka et. al— , Matsunaga et. al— , and more recently by Sergienko and 
Curnoeii. There have been a variety of microscopic proposals for the superconducting states in PrOs 4 Sbi2i^2iS£. 
It is interesting to note that none of these microscopic theory attempts to provide an origin for the two transitions 
; I ■ for cubic or tetrahedral materials. Existing microscopic theories predict the properties of the A phase but are forced 
to make ad-hoc assumptions about the origin of the transition in the B phase. This is in sharp contrast to the case 
of UPt3 where microscopic justification for the second transition exist based on weak coupling theories^. 
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In the hope of identifying a common origin to multiple transitions for heavy fermion superconductors, we apply the 
conceptual frameworks developed for hexagonal UPt3 to cubic and tetrahedral superconductors. One goal is to examine 
under what circumstances weak coupling BCS provides a microscopic description of the two transitions. We expect 
that BCS theory will capture a large contribution to the condensation energy of these heavy fermion superconductors 
and therefore provides a reasonable basis for such studies. This is partially justified by the agreement between the 
observed size of the specific heat jumps at T c and that predicted by weak coupling theories. In neither Ui_ x Th x Bci3 
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nor PrOs4St>i2 has any weak symmetry breaking fields been identified. Therefore, we begin by exploring possible 
higher symmetry groups that are weakly broken. We consider two higher symmetries groups. The first is motivated 
by weak crystal field theory of Zhitomirsky and Ueda for UPt^. In this theory, the spin-singlet superconducting 
state has an approximate spherical symmetry 5*0(3) which is weakly broken to cubic or tetrahedral symmetry. We 
next consider the possibility of weak spin-orbit coupling so that the spin-triplet superconductor has an approximate 
O x 50(3) symmetry. We find that weak coupling theory for only the latter of these two higher symmetries allow for 
one possible scenario for two transitions. The sequence of transitions found this way does not agree with experimental 
data for PrOs4Sbi2. We further show that the weak crystal field theory does allow for two transitions when strong 
coupling corrections are included. However, the condensation energy associated with the strong coupling corrections 
must be comparable to that of weak coupling theory to account for experimental results on the specific heat. These 
results indicate that a weak coupling theory does not provide an adequate description of the superconducting state 
within the context of weakly broken higher symmetry groups. 

Given the failure of weak coupling theory for PrOs4Sbi2 in the above context, we finally ask if it is possible for 
weak coupling theory to provide a description of the two transitions that agrees with the experimental results. We 
find that it is possible for such an approach to succeed and give one example of when it does. However, without the 
constraints imposed by higher symmetries, a general analysis of resulting weak coupling theories is not possible and 
requires a detailed knowledge of the quasi particles within the material. 

The paper is organized as follows. We initially provide an overview of experimental results for PrOs4Sbi2 and 
Ui-^Thjj-Beia. Then we provide an analysis that results in all the possible high temperature superconducting phases 
allowed within weak coupling theory for cubic and tetrahedral superconductors. Finally, we examine the origin of the 
second transition within the frameworks discussed above. 



II. EXPERIMENTAL PROPERTIES OF THE SUPERCONDUCTORS PROS 4 SB 12 AND Ui^TH^BEw 

The heavy fermion superconductor PrOs4Sbi2 is the first among the rare-earth filled skutterudite compounds 
showing a superconducting behavior. It undergoes a normal to superconducting transition with T c % = 1.85 K (the 
high temperature phase is known as the A phase) followed by another transition at T C 2 = 1.75 K from the A phase 
to the B phase. This second transition shows up as a pronounced anomaly in the specific heat and magnetization 
measurements 21,22 . Note that the specific heat measurements by Measson et al. have raised questions about the 
intrinsic nature of the double superconducting transition^,. Various experiments have reported a similar field 
dependance of the H C 2 curve for both the transition temperatures 2 ^. Though the symmetry and type of nodes of the 
gap structure in the A phase remains inconclusive, experiments suggest the presence of two point nodes in the B 
phase. This has been observed by the power law temperature dependance in specific heal£, thermal conductivity 2 ^ 
and penetration depth measurements^. A possible third phase transition(T C 3) has been observed as an enhancement 
of the lower critical field H c i{T) below T« 0.6 Interestingly this anomaly has not been detected by specific 

heat measurements at low temperatures. Experiments have also observed local broken time reversal symmetry in the 
superconducting phase through /iSR measurements 2 ^. 

One experimental result on PrOs4Sbi2 which has not received much theoretical attention yet provides a strong 
constraint is the low field measurement of the vortex lattice geometry. The key result is that the vortex lattice 
is not hexagonal near H c i for the field along the c axisSi. As the authors oS& have pointed out, this is a strong 
constraint because this lattice structure results from a London free energy whose form depends upon the symmetry 
of the superconduting state. The London free energy can be expanded in powers of the reciprocal lattice vector q of 
the vortex lattice. The form of the free energy is2£ 



here A is the penetration depth and m is the normalized London effective mass tensor whose components in weak 
coupling theory for a singlet superconductor is given by m~- 1 cx< VfiVfj\A(k)\ 2 >ps with vji being the i th component 
of Fermi velocity and |A(fc)| representing the gap magnitude. We are justified in keeping q terms up to second order 



near H c i because in this region q has a small magnitude being approximately given by q « « / ■¥- . If we consider that 



the gap to be invariant under a three fold rotation we find that the components of the effective mass tensor are given 
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TABLE I: The irrep and corresponding representative basis functions for even and odd parity states of cubic symmetry 



This situation would result in a hexagonal vortex lattice near H cl . No other symmetry element of the tetrahedral point 
group implies a hexagonal vortex lattice. Since the observed vortex lattice is not hexagonal near H&, we conclude 
that superconducting state does not contain a three fold symmetry axis. 

The alloy Ui-^Th^Beia which has cubic point group symmetry shows two second order superconducting transition 
for concentration of thorium exceeding x = 0.018 in specific heat measurements^. A pronounced peak has been 
observed in the ultrasonic attenuation for longitudinal sound propag ated along a [100] direction at T = T c ^. Mea- 
surements by Bishop et al. found a similar behavior for longitudinal sound propagated along the [111] direction- 2 . 
Significant anomalies have been observed in H c i (T) at T — T c ^^ and the zero field ^SR line width shows a marked 
increase as T decreases below T C 2 in samples with x rj 0.032 3 ! Finally Lambert et al. found differences in the pressure 
dependance of T c for x < 0.018 and x > 0.018^. 

III. THEORY OF THE HIGH TEMPERATURE SUPERCONDUCTING PHASE 

For all the theories we consider, the initial transition into the superconducting state is characterized by either cubic 
or tetrahedral symmetry. Therefore, this transition is described by a single order parameter symmetry. Here, we 
provide a general analysis of the possible superconducting states that weak-coupling theory allows for this phase. A 
similar analysis of this problem has appeared recently by Kuznetsova and Barzykin 36 . However, we find that weak 
coupling theory provides even stronger constraints than found in this work. As a consequence, some of the phases 
found by these authors are ruled out. 

The cubic group has ten irreducible representations (irreps) that are listed in Tabled The tetrahedral group has 
one three dimensional irrep and three one dimensional irreps, two of which are complex conjugate and combine to 
form a single irrep when time reversal symmetry is present. In this Section we study the superconducting phases of 
materials with cubic and tetrahedral point group symmetry whose basis functions transform as a multidimensional 
representation. 



The superconducting gap is given by A(fc) = [iJj(k)ao + d(k) ■ cf\ia v where ip(k) — ip(—k) is the even-parity spin- 
singlet component, d(k) = —d(—k) is the odd-parity spin-triplet component, and <x = (a x , a v , a z ) represent the Pauli 
spin matrices. Among the various irreps of the point group symmetry there is one, say T, which gives the highest 
transition temperature. The superconducting state can be written as a linear combination of the basis functions of 
this representation A(r, m; k) 

A(fe) =^77(r,m)A(r, m; k) 

m 

Here r](T, m) are in general complex and act as the order parameter. For a single representation the fourth order free 
energy for the superconducting state can be written as an expansion in ry(r, m) 

F r (T, V ) = F (T) + a £ ^(r, m)| 2 + f r [ V (T, m) 4 ] 

m 

with a — cto[T/T e (T) — 1], and the fourth order energy contains all terms which are invariant under G x R x U(l). 
The normal state of the system is represented by the free energy Fq (T) . Here G is the crystal point group symmetry, 
R and U(l) are the time reversal and gauge symmetry groups respectively. The fourth order terms in the Ginzburg 
Landau theory are characterized by several coefficients (3- L which are arbitrary in a general theory but are constrained 
in the weak coupling limit. 
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For the spin triplet systems the fourth order free energy within weak-coupling theory can be evaluated from the 
average 

(|d(fe)| 4 + |q(fe)| 2 ) 

with q = id(k) x d(k)* . The q vector is zero for unitary states and takes finite values for non unitary states. Since 
|<7| 2 gives a positive fourth order contribution to the free energy for non unitary states it is unusual for these states 
to be preferred stable states within a weak coupling theory. For the spin-triplet representations, sixth order terms 
in the free energy are required to remove a residual degeneracy and completely specify the solution. The sixth order 
free energy is given by 

i(|d(fe)| 6 + 3| 9 (fc)| 2 |d(fc)| 2 ). 
For spin singlet systems the fourth order free energy will be given by 

(\m\ 4 )- 

We now analyze the possible superconducting ground states for cubic symmetry within the weak coupling limit. 



A. Two Dimensional Representation 



The cubic group contains one two dimensional representation for even parity states with the gap function 

4>(k) = mfiE 3 (k) +m.f2E g {k) 

where fiE g (k) and f2E g (k) form a basis for the E g irrep. For the odd parity states the representative basis functions 
f 1E and f 2E are given in Tableland the gap function is given by 

d(k) = rnf XEu {k)+r) 2 f 2Eu {k) 

The general free energy can be expressed as an expansion in the order parameter 771 and 772 as 

T 



F = a 



1 

T c 



(M 2 + \m\ 2 ) + yIII^oA^M 2 + M 2 ) 2 + - nlmf] (2) 



Where «o = 1 m the weak coupling limit. The weak coupling values of the fourth order coefficient for the even parity 
states are 

A=3/3 2 =</ 1 4 £g (fc)> 

with the bracket meaning the average over the Fermi surface. All averages in the current and future discussion are in 
a normalized form such that < f 2 {k) >= 1 where f(k) represents the basis function. If we minimize this free energy 
we find that since f3i > 0, the phase uj 2 (l,i) will minimize the free energy. This is true for arbitrary Fermi surfaces 
and gap basis functions. This phase belongs to the superconducting class 0(Z?2)^ and has point nodes along the 
cube diagonals. The gauge factor u> 2 has been multiplied to keep the notation consistent with Ref. 36. 
For the odd parity irreps the weak coupling values for the fourth order coefficients are 

Pi = 3 < / 2 (fc)/ 2 (fc) > (x + l),/? 2 =< f 2 (k)f 2 (k) > (x - 7) 

where x =< f^{k) > / < / 2 (fc)/ 2 (fc) >, and [f x (k), f y (k), f z (k)] form an arbitrary basis of T\ u symmetry and have 
the same rotation properties as the vector k. If we minimize this free energy within the weak coupling limit, the non 
magnetic phase v = (1, 0) belonging to the superconducting class D4 <g) R is found to be stable for x < 7 (note that a 
residual continuous degeneracy remains for which any real combination of the two components minimizes the fourth 
order free energy, this is lifted by sixth-order terms). For x > 7 the magnetic phase v — u> 2 (l, i) is stable. Again the 
gauge factor uj 2 has been multiplied to keep the notation consistent with Ref. 36. This phase belongs to the 0(Z?2) 
superconducting class and contains point nodes along the cube diagonals. It is interesting to find a non unitary phase 
that is stabilized within weak coupling theory in zero applied field. The second non-magnetic phase v — (0, 1) appears 
to be prohibited in a weak coupling theory by the sixth order term in the free energy 

-i,\m\ 2 \^l-vl\ 2 



5 



where 73 is given by the weak coupling value 



73 = ^ < f!(k) + 2/J(fc)/ a 2 (fe)/ 2 2 (fc) - 3/ x 4 (fe)/ y 2 (fc) > 



we numerically find 73 > for a variety of basis functions and Fermi surface structures, but we could not prove this 
analytically. 

Tctrahedral symmetry does not change the structure of the w 2 (l,i), but does introduce an additional sixth order 
term in the free energy that modifies the v = (1, 0) phase. This term is 

^ - tt(u\ &lu\l -f2/r,\ f2(r,\\ ^ fo /~, „* 1 „ |4 , \„ |4 ni„ |2i |2\ /-„3,„*3 . ^3^*3\ 
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< m)fy(k)(m) #(*)) > [3(mv*2 + mv*i)(\m\ 4 + \mV - ^H^r) - (vM* + vMl > 



as a result the stable A phase ground state is given by (fa, (fa), where both fa and fa are real (note that there is no 
continuous degeneracy in this phase 

B. Three Dimensional Representation 

The Free energy for the three dimensional representation can be written as 
F = a 

+Mpi+p 2 2+pi)(pr+pr+pr)+M\pi\' + h 4 + i P3 i 4 )] 

(3) 



T 



(M 2 + Iwl 2 + \m\ 2 ) + 16 7 ff' )2 JVoA 4 [ft(| P1 | 2 + | B | 2 + | B | 2 ) 2 



The weak coupling values of the coefficients are 

ft = 2/3 2 = 2 < f? T2g (k)fi T2g (k) >,/% =< ff T2g (k) > -3 < f? T2g (k)fl T2g (k) > 

The functions [/it 2 {k),f2T 2 (k),f3T 2 (k)] form a basis for the T 2s irrep. If we define the parameter space in terms 
of one free parameter 

</m„(fe)> 



<ftrJk)fi T Jk)> 

we find that for x < 3, the state 0) is the ground state. This state belongs to the class D^E) and contains 
line nodes in the z = plane. For x > 3, the phase (1,ui,uj 2 ) is stable high temperature phase. This state belongs 
to the class D^{E) and contains point nodes. The boundary x = 3 corresponds to a spherical Fermi surface. The 
weak coupling values provide a tighter restriction on the allowed phases than found by Kuznetsova et. al. and result 
in only magnetic states being the stable weak coupling phases. In particular, it should be noted that x > 1 for any 
choice of basis functions, and this constraint rules out any non-magnetic phases. 
For the odd parity states the weak coupling values are given by 

[h = X + 3, fa = 7: , P3 = n 



where again a; =< f*(k) > / < f%(k)fy(k) >. 

We find that the state (1, 1, 1) is stable for x > 3 whereas (1, 0, 0) is stable for x < 3 for both T\ u and T 2u irreps. The 
boundary x = 3 corresponds to the spherical Fermi surface. 

All the possible solutions for the high temperature phase within the weak coupling limit are listed in table |n| For 
the tetrahedral symmetry the analysis will be similar to the cubic at the fourth order as tetrahedral and cubic group 
have the same invariants. The difference in the tetrahedral and cubic invariants at the sixth order results in the state 
( 1 771 1 , * 1 772 1,0) belonging to D2(E) symmetry group giving the ground state for the tetrahedral symmetry rather than 
the (l,z,0) state belonging to the D±{E) symmetry group giving the ground state for cubic symmetrjiii. 

In view of the controversy between the extrinsic versus intrinsic nature of the transition in Pr04Sbi2, we briefly 
consider that there is single superconducting transition and ask if any of the above stable phases can explain the 
experimental properties observed at low temperatures. In particular, the ground states listed in table ITT1 which may 
explain the observed properties are the (1,0) state of E u , (1, i, 0) for T2 g , and the (1,0,0) state of T\ u and Tm- 



6 



Representation 


State 


Nodes 


Symmetry Class 


Eg 




P 


O(Ai) 


E u 


(1,0) 




D4 g> 7? 


Fju 




JD 
I 


U(L>2 ) 


T 2g 


-~(1,*,0) 


L 


D A (E) 




75(1, w,w ) 


P 


£> 3 (#) 


Tiu 




P 






(1,0,0) 


P 


D 4 {D 2 )®R 




^a.M) 


P 


D 3 (C 3 )®R 




(1,0,0) 


P 


D 4 (C 4 ) ®R 



TABLE II: Stable high temperature phases in the weak coupling limit. The corresponding representative basis functions are 
listed in table [I] Here u = exp[27ri/3] 



These states have a gap structure with point nodes or are highly anisotropic. In addition, each of these states 
break tetrahedral symmetry which would result in a distorted vortex lattice structure. Note that the observation of 
an increased muon relaxation rate in the superconducting phases^ does not imply that time reversal symmetry is 
globally broken, but only locally broken—. We take this to imply that the order parameter must be multi-component 
but that ordered phase need not break time reversal symmetry globally. 



IV. WEAKLY BROKEN SO(3) THEORY FOR SPIN SINGLET STATES 



We now turn to the role that higher symmetries may play in giving rise to two superconducting phase transitions. 
In this section, we consider the weak crystal field theory for which there is an approximate SO(3) for spin-singlet 
superconductors. Such an approach was proposed to explain the phase diagram of UPt3 by Zhitomirsky, and Ueda 8 . 

We find the possible superconducting transitions for a state in which the spin singlet Cooper pairs are in the I = 2 
channel. Due to the effect of a weak crystal field, the five- fold degenerate I — 2 irrep of 50(3) split into E g © T 2g of 
the cubic group. The free energy for the I = 2 irrep of SO(3) has been found by Mermin and Stared, incorporating 
the weak cubic field gives 

/ = ai (\ m \ 2 + \ m \ 2 ) + a 2 (\ Pl \ 2 + \ P2 \ 2 + \ P3 \ 2 ) + [3i\TrB 2 \ 2 + l3 2 (TrB*B) 2 + f3 3 Tr(B 2 B* 2 ) (4) 

here B is a 3 x 3 traceless symmetric complex matrix given by the I = 2 order parameter ip(k) = ^2^ u B^k^k^. 
Note that since the symmetry breaking is weak, spherical symmetry is broken by the second order term only. The 
weak coupling limit corresponds to the special case^ /3 2 = 2/3%, (3 3 = 0. The magnitude of (3 3 is only due to strong 
coupling effects and is of order T c /Ep ■ The general form of the gap is 



^(fc) = ^( 2fc 2 2 ~ k l ~ k l) + ^=( fc ' ~ k 2 y ) + V2 Pl k x k y + V2 P2 k y k z + V2 P3 k z k x 

here r\ — (771,772) transform like the E g representation and P = { P i,p 2 , P3 ) transform like the T 2g representation of 
the cubic group. The components of matrix B can be written as 

u ^1 , d r /i d 2r li d n n 



We now look for the various transitions into the B phase which can result in further stable superconducting transitions. 
The form of this theory is simple enough to enable us to perform a general analysis beyond the weak coupling limit. 
We will therefore look for all possible transitions from all possible stable A phase solutions listed in—. From the values 
of the fourth order energies in Table ITTT1 we find that for (3 2 > 0, f3 3 > the states P = (1, i, 0) and P = (1, oj, cu 2 ) will 
tend to stabilize deep in the B phase when the second order coefficients can be ignored whereas for (3 2 > 0, j3 3 < 
the state rj = u> 2 (l, i) tends to stabilize. Within the weak coupling theory these states do not give any stable second 
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TABLE III: The fourth order free energies {Fa) for possible A phase representation (r) for spherical symmetry. Here r\ — (771 , 772) , 
P = (pi,P2,p-i) transform as irreps of E g and T2 g respectively of cubic symmetry. Here u> = exp[2-7ri/3] 



order transitions into the B phase. If we include strong coupling effects (ft ^ 0) we find that there is only one second 
order transition for ft > 0, ft < from the state P = -^(1, i, 0) in the A phase. This transition corresponds to a B 

phase given by a linear combination of P = ^(1, i, 0) and r\ = j^e^l, 0) with a transition temperature 

4ft 

T cB = T cA + 3(1 + -^)(T cA - T<) 

ft 

where T < is the transition temperature corresponding to pure r\ = (771,^2) state. This state is highly anisotropic and 
gives a distorted vortex lattice structure. It is therefore a possible transition sequence for PrOs4Sbi 2 . The specific 
heat jump ratio between the transition at the B phase to the transition at the A phase is 

C B ^T cB ft 2 (ft+ft/4) 

C A T cA 144(1 + ft + ft/2)(ft + ft/3) 2 

This transition gives a specific heat jump ratio of the order ft 2 which is negligible close to the weak coupling limit. It 
is also interesting to note that this transition corresponds to a change in penetration depth of order ft which can be 
a significant change to observe in an experiment. A similar explanation may hold for the T c3 in PrOs4Sbi2 observed 
in penetration depth measurements below T=0.6K but not yet observed in specific heat measurements. 



V. WEAKLY BROKEN O x SO (3) THEORY FOR SPIN-TRIPLET STATES 

We will now analyze the transition for the spin triplet states within a weak coupling theory. We consider the effects 
of a weak spin orbit coupling and also include the crystal field with cubic symmetry and allow the spin channel to 
be isotropic. The irreps of the symmetry group are given by the combined group Oh x 50(3) x R. If we consider a 
weak spin orbit coupling in our system, the basis functions split up into four different irreps. We write the vector gap 
equation in terms of these irreps as d(k) = Sid ri (fc) where the components are given by 



,{k) -yf y (k)} +V2^[-2zf z (k)+xf x (k)+yf v (k)} 

(k) - zf y (k)\ + P2 l=[xf y (k) - yf x (k)\ +P3 ±[zf x (k) - xf z (k)} 

(k) + zf y (k)} + q 2 ^=[xf v (k) + yf x (k)] + g 3 -L[z/ x (fe) + xf z (k)] 

(5) 

the components of the vector gap equation can be written in terms of the basis d{k) = ^ . rjjfj(k), and we will use 
both bases for convenience. 

For a weak spin orbit coupling, we get a free energy of the form 

F = «i|A| 2 + a2 (M 2 + M 2 ) + a 3 (N 2 + \P2\ 2 + \ P3 \ 2 ) + MM 2 + M 2 + M 2 ) + ftdryj 4 + \v y \ 4 + \v z \ 4 ) 
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TABLE IV: The fourth order free energies(J4) for possible A phase irreps(F) and their corresponding weak coupling val- 
ues(W.C). The weak coupling values are in units of < fx(k)fy{k) >. The value of x is given by the ratio — f^u^J'L — ■ Note 
that x=3 corresponds to spherical symmetry, uj — exp[2Tri/3]. 



+fo(\vl\ 2 + W y ? + \vT) + k(\v x -V v \ 2 + \riy-ri z \ 2 + \r, x .r, z \ 2 ) + ft [(t^;) 2 + (r^.r,;) 2 + {r]y.rf z f + c.c] 
+fc(\ri x .r 1 * y \ 2 + \VyVl\ 2 + \V X -Vl\ 2 ) + fo[(n x ?{rfy? + (^(ril) 2 + (v y ) 2 (vl) 2 + «] 

+Pr(\ri x \ 2 \r 1y \ 2 + \ Vy \ 2 \rif + \r ]z \ 2 \r 1x \ 2 ) (6) 
The weak coupling values of the normalized coefficients are 

ft = -2ft =< g(k) >,ft = \fo = -2ft - ift = -ift =< f 2 x {k)f 2 y {k) > 



It may be seen from Table ITTT1 that for x > 3 the states P = (1, 1, 1) and Q = (1,1,1) will minimize the fourth order 
terms whereas for x < 3, v — (1,0), P = (1,0, 0),Q = (1,0,0) and equivalent states minimize the fourth order free 
energy. At x = 3 we have a spherical fermi surface. 

To understand if phase transitions are possible we compare the states that minimize the second order term with 
those that minimize the fourth order terms. If these arc different, then a transition is possible. However many of 
these transitions are first order. For example the non unitary state v — aj 2 (l,i) which is stable in the A phase for 
x > 7 undergoes a first order transition to a mixed state with v = ui 2 (l,i) and P = — u>, lo 2 ). 
We find that there is only one stable second order transition into the B phase. This instability is from P — p(l, 1, 1) 
to the combination of P — p(l,l,l) and Q = g(l,l,l), where p and q are the order parameter values. Within 
the assumption of a small spin orbit coupling we write the general form of the gap for the states P = (1,1,1) and 
Q = (1,1,1) as 

A(k) = cos^^ftfe)^^^)]^^^)^^^)]^^^)^^^^^)]) 

sm6([xf y (k) - yf x {k)\ + [zf x (k) - xf z (k)} + [yf z (k) - zf y (k)]) (7) 

Here 6 acts as the order parameter. If we assume that the P = (1, 1, 1) state has the highest transition temperature, 
there is a transition from a 9 = state in A phase for which 9 becomes non zero and grows towards a fully gapped 
system at 9 — 7r/4 (at which the gap has the form xf y (k) + yf z (k) + zf x (k)). This gives a second order transition 
temperature 

TcB = T cA - - (x + 3)(T cA - T<) 

and the ratio of jumps in the specific heat at the transition temperatures of the A and the B phases is 

cb _ TcB 
c A T c a{x + 3) 

It should be pointed out that for tetrahedral group there will be an additional bilinear coupling term in the free energy 

a m (P* Q + P Q*) 
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r 


F 4 


r? = (l,0),P = 


ft 


2 


2ft/3 


= o,p = ^(i,»,o) 


(ft +ft)/2 


V = 0,P = (1,0,0) 


ft 


*» = 0,P = ^(l,w,w 2 ) 


ft/3 + ft 


i7=^(0,l),P = ^(0,0,<) 


l/4(ft + ft + 2ft -4ft) 


rj = 0,P= ^(1,1,1) 


ft/3 + 2ft 



TABLE V: The fourth order free energies(_p4) for irreps of stable states(F) in the A phase 

which would smear out the transition. 

Though we get two transitions in this case, the sequence does not satisfy the observed physical properties in 
the skutterudite PrOs4Sbi2- Owing to the tendency of this system towards a fully gapped state such a transition 
would give a vortex lattice structure which becomes hexagonal at the lower temperatures which is contrary to the 
experimental observations where the distortions from a hexagonal structure increase at low temperatures^. In 
addition this state does not satisfy the nodal structure of the gap in the B phase as observed in thermal conductivity 
and magnetization measurements. 



VI. WEAK COUPLING ACCIDENTAL DEGENERACY THEORIES 

Here we consider the accidental degeneracy between the T 2g and E g irreps as an example, since as explained at the 
end of this section there are many theories that have multiple transitions. 

In weak coupling theory we can write the fourth order free energy by evaluating the average 

< \lP(k)\ 4 > = < \rnf 1Eg (k) +mf2E g (k) +PlflT 2g (k) +P2f2T 2g (k)+p 3 f 3T2g (k)\ 4 > 

In the above expression r/ = (771,772) and P — {pi,p 2 ,p 3 ) are the order parameter values and the representative basis 
functions / are assumed real. The free energy expression is 

F = ai (\ Vl \ 2 + \ m \ 2 ) + a 2 (\ Pl \ 2 + \ P2 \ 2 + \ P3 \ 2 ) + Mhil 2 + M 2 ) 2 + \(viv; - V2V*i) 2 } 

+ft(M 4 + b 2 | 4 + b 3 | 4 ) 

+/3 3 [4(|pi| 2 | P2 | 2 + \P2\ 2 \ P3 \ 2 + bslVl 2 ) +p\pf + P%pf +Ph? + cc] 

+ft/2[4(M 2 + M 2 )(N 2 + b 2 | 2 + bsl 2 ) + ((vi + Vi)(p*i 2 +P*2 2 +Pf) + c.c)] 

-/3 5 /2[4(M 2 - |r72| 2 )(2|p 3 | 2 - | P i| 2 - |p 2 | 2 ) - + vlm){\pi\ 2 - |p 2 | 2 ) + 

[{ill ~ r&){2pf-p?-p?) - 2V3 m r l2 (pf - p* 2 ) + c.c}} (8) 

Within a weak coupling theory the normalized coefficients are given by the following cubic averages 

0i =< f? Bg (k) >,fo =< A 4 T2s (fc) =< f 2 T 2g (k)f 2 T2g (k) >,ft =< [f 2 1Eg (k) + fl Eg (k)}f 2 T Jk) > 

f% =< [fiE 9 (k) - fl Eg (k)}f 2 T Jk) > 
The spherical fermi surface corresponds to the special case 

ft = /? 2 - 3/3 3 - |/3 4 , /3 5 = 

From Table [V] we find that for ft < the lowest 3 fourth order energies in increasing order correspond to states 
P = 0), P = (1, w, to 2 ), r/ — (1, i). We do not find any second order weak coupling transitions within this range. 
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FIG. 1: (color online) Gap structure in the B phase for the P = \r\ (1, i, 0) + r\ = e ,7r ' 4 (|?7i |, I772 1 e l7ry/2 ) phase. 

For ft > the lowest 3 fourth order energies in increasing order correspond to states r\ = (l,i),P = (1,cj,cj 2 ),P = 
0). In this case we find only one stable second order transition from the state P = 0) in the A phase to 

a B phase where it mixes to the state rj = e l7r / 4 (|r7i|, \r]2\e l7T ^ 2 ) with a transition temperature 

T c b = T c a — — \T C A ~ T<) 

1 2/3j-V7p 5 y <J 

The specific heat jump ratio at the transition temperatures is given by 

cs = Tcs ( 7(ft 2ft + V7f3 5 ) 2 \ 

ca T cA y 12ft(ft + ft) - 7(2ft - v^ft) / 

We find that a specific heat jump ratio that is comparable to observed value for the skutterudite PrOs4Sbi2 which 
is about oneSican be obtained if the gap functions contain substantial cubic anisotropy. In addition this situation 
will result in a anisotropic state with two fold degeneracy. This sequence of phase transition would also result in a 
distorted vortex lattice structure owing to the large anisotropy of this state as shown in Figure 1. 

We find that in many cases, the A to B transition is first order. The reason for this is as follows: for two order 
parameters ip and r\ the free energy takes the form 

F = ai|Vf + a 2 \v\ 2 + PM 4 + foM 4 + PmiWM* + MVV 2 + »?V 2 ) 
Let (ip, rj) = (\tp\-, |??|e z ^) and minimize with respect to <fi. The the free energy becomes 

F = aM\ 2 + a 2 M 2 + ft + f3 2 \r,\ A + /3 m |r?| 2 |V| 2 

where f3 m = ft m i — 2/3„ l2 ~< / 2 (fc)/, 2 (fc) >. Here we have assumed the basis functions fy(k) and f v (k) to be real. 
If there is a second transition then it is second order transition when (3^ < Aflifa, otherwise it is first order. In our 
calculations /3 m is too large, which leads to the first order transitions. This is a consequence of the functions f^p(k) 
and f v (k) that have been chosen. However we can get a second order transition by considering a two band theory for 
which f,j,(k) is large and f v (k) is small in one band while f v (k) is large and f^(k) is small on the other. Then the 
condition 0^ < 4/3i/3 2 will be easily satisfied and a second order A->B phase transition can exist for almost any two 
different order parameter irreps. 

VII. CONCLUSIONS 

We have considered microscopic theories of unconventional superconductivity in cubic and tetrahedral superconduc- 
tors. We have identified the stable weak-coupling unconventional superconducting states that belong to a single irre- 
ducible representation and have highlighted which of these can describe the low temperature properties of PrOs4Sbi 2 . 



11 



We have further examined theories for two intrinsic superconducting transitions in PrOs4St>i2. We have found that 
a theory for which the two transitions are due to a weakly broken SO(3) symmetry for spin singlet Cooper pairs 
cannot give rise to two transitions in the weak-coupling limit. However, it is possible for such a theory to produce 
two transitions that agree with the experimental properties of PrOs4Sbi 2 only when extended to the strong coupling 
limit. We further find that for spin-triplet Cooper pairs, weak spin-orbit coupling in the weak-coupling limit does not 
give rise to two superconducting transitions that agree with the experimental properties of PrOs4Sbi2. Finally, we 
consider an example of a weak coupling theory that does not assume an approximate higher symmetry, but is based 
on the accidental closeness of the transition temperatures for two different representations. This example is able to 
describe the observed properties of PrOs4Sbi2. 

This work was supported by the National Science Foundation grant No. DMR-0381665. 
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